Bonn-TH-2012-033 



The Stability Of R-Parity In Supersymmetric Models Extended 

By U{1)b-l 

J. E. Camargo-Molina0 B. O'Leary0 and W. PorocQ 

Institut filr Theoretische Physik und Astronomie, Universitdt Wiirzhurg 
Am Huhland, 97074 Wuerzhurg 

F. StautI!] 

Bethe Center for Theoretical Physics & Physikalisches Institut der Universitdt Bonn, 

53115 Bonn, Germany 

Abstract 

We perform a study of the stability of i?-parity-conserving vacua of a constrained version of the 
minimal supersymmetric model with a gauged U{1)b-l which can conserve i?-parity, using homo- 
topy continuation to find all the extrema of the tree-level potential, for which we also calculated the 
one-loop corrections. While we find that a majority of the points in the parameter space preserve 
i?-parity, we find that a significant portion of points which naively have phenomenologically ac- 
ceptable vacua which conserve i?-parity actually have deeper vacua which break i?-parity through 
sneutrino VEVs. We investigate under what conditions the deeper i?-parity-violating vacua ap- 
pear. We find that while previous exploratory work was broadly correct in some of its qualitative 
conclusions, we disagree in detail. 
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I. INTRODUCTION 



With the discovery of a resonance at 126 GeV at the LHC [HE], the interpretation that 
it is the Higgs boson of the standard model (SM) raises the issue of the stabihty of the SM 
vacuum [3]-[7j. Since there is only one field in the SM that can possibly have a non-zero 
vacuum expectation value (VEV) (assuming Lorentz invariance), finding the minima of the 
potential energy is straightforward, though of course evaluating it to the accuracy required 
is quite involved [5-7J. 

Many extensions of the SM introduce extra scalar fields, and thus could in principle have 
many more VEVs. In particular, the minimal supersymmetric extension of the SM (MSSM) 
not only introduces complex scalar partners for each SM fermion, it also introduces a second 
Higgs SU{2)l doublet. Even only allowing for the neutral components of the Higgs doublets 
to gain non-zero VEVs, the vacuum structure is non-trivial. Fortunately, as a particular 
case of a two-Higgs-doublet model, there is a proof that at tree level, there are at most two 
minima that conserve electric charge [8j , and there are analytic formulae to determine which 
vacuum has a lower energy at tree level [9j . 

However, the possibility of VEVs for the scalar partners of the SM fermions, the sfermions, 
is a troubling prospect, as all but the sneutrinos are charged under SU{?>)c and/or [/(l)e^, 
and thus any VEVs for these scalars are ruled out by these being good symmetries of the 
vacuum. Unfortunately, the potential as a function of several scalar fields is so complicated 
that even at tree level, finding the global minimum is highly non-trivial, in general [10]- 
[13]. The only tractable approaches thus far have been to investigate the "most dangerous" 
possibilities of stop or stau VEVs [H [T4l [T5] . 

Sneutrino VEVs, on the other hand, are not necessarily a problen|^ Of course, they would 
break i?-parity, leading to mixings between SM and SUSY particles carrying the same spin 
and SU{?))c x U{l)em quantum numbers, see e.g. [17j and refs. therein. However, in the 
absence of baryon-number-violating trilinear couplings in the superpotential, proton decay 
would still be perturbatively forbidden. If i?-parity is conserved, the lightest supersymmetric 
partner (LSP), if uncharged under SU (3)c and U (l)em5 could be a perturbatively stable dark 
matter candidate pSl HF]- If i?-parity is not conserved, the LSP is not absolutely stable, 
^ For example, in the case of bilinear i?-parity, the sneutrino VEVs have to be small to explain neutrino 
data implying that in this model essentially the conclusions concerning the minima of the potential can 
be taken over from the MSSM 
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but, depending on the size of the sneutrino VEVs, a smaU amount of i?-parity violation 
stiU aUows for a metastable LSP to make up the observed dark matter. However, even if 
the neutrahno is not sufficiently long-lived, another particle like the gravitino [20l - [24] or the 
axion [241426] could be long-lived enough to serve as dark matter. 

The conservation of baryon number B and of lepton number L is only accidental in the 
SM, and this conservation through the imposition of i?-parity may appear to be rather ad 
hoc in the MSSM, since the only difference between a lepton SU{2)l doublet and the Higgs 
SU{2)l doublet responsible for mass for the down-type quarks Ha is that is assigned 
L = instead of 1. StiU, noting that i?-parity is defined as = (_i)(3(b-l)-2.) 

if lepton number is integer or half-integer for every field, one may consider a symmetry based 
on B — L as a motivated extension. 

Adding as an extension to the SM [27H3T] or the MSSM [32ti34] has been con- 

sidered quite a lot in the literature, as in addition to motivating i?-parity, it leads to a 
rich phenomenology, and, in some models, naturally leads to type I see-saw mechanisms for 
neutrino masses. 

The most minimal of such extensions of the MSSM requires superfields with right-handed 
neutrinos as fermionic components to cancel gauge anomalies, and their scalar components 
provide the VEVs necessary to break U{1)b-l as required by phenomenology. Already such 
a model has eight complex scalar fields that in principle could have non-zero VEVs without 
breaking SU{?>)c or [/(l)em, with many possible minima. 

Conserved i?-parity is still attractive from a phenomenological point of view, since there 
are tight limits on many i?-parity- violating parameters [351437] . One can extend the neutral 
sector with further superfields with VEVs which break U{1)b-l without breaking i?-parity 
[321 l3lj . However, given the complicated structure of the potential, it is prudent to be 
concerned about whether such extensions really do conserve i?-parity at their global minima 
over relevant parameter ranges. 

In the following sections, we investigate how stable i?-parity-conserving vacua are in such 



a minimal [/(l)^_L-gauged extension of the MSSM. In sec. II, we define the model. In sec. 
Tin we describe how we approach the difficult problem of minimizing such a complicated 
potential, including going to the full one-loop potential, since it is well-known that there can 



be rather large corrections in the Higgs sector of the MSSM [38146T] . In sec. IV, we present 



the results of two scans over phenomenologically interesting regions of the parameter space. 



and in sec. |IV A| we compare our results to previous investigations in the literature. Finally, 
we sum up in sec. [V, 



II. THE MODEL 

A. Particle content and superpotential 

There are several ways to extend the MSSM by U{1)b-l- Here we choose the minimal 
extension which allows for a spontaneously broken U{1)b-l without necessarily breaking 
i?-parity. This requires the addition of two SM gauge-singlet chiral superfields carrying 
B — L which may develop VEVs, as well as the addition of three generations of superfields 
containing right-handed neutrinos. We refer to this model as the BLSSM. 

The BLSSM appears to be a relatively straightforward extension of the MSSM, yet it 
has a rich phenomenology, with a Z' boson (with prospects for the LHC discussed in j62]). 
Majorana neutrinos with see-saw masses, several qualitatively new dark matter candidates 
with respect to the MSSM ^J, and a rich Higgs boson sector [64j. It also has many 
interesting technicalities [65] , which are fully explained in [66j . 

However, all the above works assume that i?-parity is conserved. Indeed, the model 
was constructed as the minimal extension of the MSSM that includes U{1)b-l as an extra 
gauge symmetry, while allowing it to break spontaneously without necessarily violating the 
conservation of i?-parity. As mentioned in the introduction, R-sneutrino VEVs could break 
U{1)b-l at the cost of losing i?-parity, and aspects of the LHC phenomenology in the case 
of broken i?-parity are discussed in [671 EH] • The purpose of this work is to investigate how 
robustly i?-parity is conserved for parameters of phenomenological interest. 

While we refer the reader to Ref. [66] for a full description of the BLSSM and its imple- 
mentation in SARAH [691472] and SPheno [73l[71], for reference we recap the gauge symmetries 
and particle content of the model here. 

The model consists of three generations of matter particles including right-handed neu- 
trinos which can, for example, be embedded in S0{1Q) 16-plets. For convenience, we refer to 
their scalar components as R-sneutrinos. Moreover, below the GUT scale the usual MSSM 
Higgs doublets are present, as well as two fields r] and f) responsible for the breaking of 
the U{1)b-L' Furthermore, vj is responsible for generating a Majorana mass term for the 
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TABLE 1. Chiral superfields and their quantum numbers. 



right-handed neutrinos and thus we interpret the B — L charge of this field as its lepton 
number, and hkewise for 7), and caU these fields bileptons since they carry twice the lepton 
number of (anti-)neutrinos. We summarize the quantum numbers of the chiral superfields 
with respect to the model's gauge group U{1)y x SU{2)l x SU{?>)c x U{1)b-l in Tab. [l} 
The superpotential is given by 



+ Yi^ vILj Hu-^i'vfi+ Y'J vlf^v] 
and we have the additional soft SUSY-breaking terms: 



(1) 



2|^|2 



^SB =^MSSM - ^B^B'^BB' " T^^B'^B'^B' " rn\\vif - m"^\r\ 



1 

2' 



(2) 



are generation indices. Without loss of generality one can take and B^i to be real. 
The extended gauge group breaks to S\J{^)c ® f/(l)em as the Higgs fields and bileptons 
receive vacuum expectation values (VEVs): 



H^=-^ {vd + cTrf + i(j)d) , 
?7 = (v^ + (T^ + i(/)^) , 



V = 



V2 



{vu + <yu + i4>u) 



V2 



(3) 
(4) 
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We define tan^S' = Vj^/vf^ in analogy to the ratio of the MSSM VEVs (tan^S = Vu/va). For 
certain parameter combinations a spontaneous breakdown of i?-parity can occur as also some 
or all sneutrinos can obtain non- vanishing VEVs. We denote the VEVs for the sneutrinos of 
the SU{2)l doublets Li by Vl^i and those of the SU{2)l singlet sneutrinos by vr^i^ with 
i = 1,2,3. 



B. Constrained model 



We will consider in the following a scenario motivated by minimal supergravity assuming 
a GUT unification of all soft SUSY-breaking scalar mass parameters as well as a unification 
of all gaugino mass parameters 

2 2 2 2 2 /r\ 

2-1 2 2 2 2 2 2 /r\ 

ttIqI =m^ = rrijj = ttiq = = = m~c (oj 
Ml/2 =Mi = M2 = Ms = M^, (7) 

Also, for the trilinear soft SUSY-breaking coupling, the ordinary mSUGRA-inspired condi- 
tions are assumed 

Ti = AoYi, i = e,d,u,x,iy . (8) 

Furthermore, we assume that there are no off-diagonal gauge couplings or off-diagonal 
gaugino mass parameters present at the GUT scale, motivated by the possibility that the 
two Abelian groups are a remnant of a larger product group which gets broken at the GUT 
scale, as discussed in [66j, to which we refer the reader for details of the gauge coupling 
structure. 

In addition, we consider the mass of the and tan as inputs and use the following set 
of free parameters 

mo, Ml/2, ^0, tan/3, tan/3', sign(/x), sign(/x'), m^/, and 1^,. (9) 

Yjy is constrained by neutrino data and must therefore be very small in comparison to 
the other couplings in this model, as required by the embedded TeV-scale type-I seesaw 
mechanism, so we take Y^^ = 10~^5'^^ as its precise structure does not affect any of our 
conclusions. Y^ can always be taken diagonal and thus effectively we have 9 free parameters 
and 2 signs. We denote this constrained version of the BLSSM as the CBLSSM. 
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III. MINIMIZING THE POTENTIAL 



The primary question that we wish to answer is: within the CBLSSM, what portion 
of parameter points which have a phenomenologicaUy acceptable, i?-parity-conserving local 
minimum of the potential turn out to have a global minimum with different VEVs, and are 
there any patterns in the parameters of such points? By "phenomenologicaUy acceptable" , 
we restrict ourselves to parameter points which have a local minimum with the expecta- 
tion values of the Higgs doublets leading to the correct values for rriw and m^, there are 
no tachyons, and no charged or colored scalar has a VEV. Hence any minima that are 
''phenomenologicaUy acceptable" and i?-parity-conserving have expectation values for the 
Higgs doublet fields Hu and the bilepton fields 77,7), and no other field has a non-zero 
expectation value. 

Our method to address this question was as follows: we performed a random scan over a 
range of input parameters, constrained to having phenomenologicaUy acceptable, i?-parity- 
conserving local minima at expectation values for the Higgs doublet and bilepton fields 
given as input; we then found all the extrema of the tree-level potential for each parameter 
point allowing non-zero sneutrino VEVs; we evaluated the loop-corrected potential at these 
points (ensuring that the corrections did not move the minima significantly) and selected 
the deepest, and thus can say if the minimum chosen as input was not the global minimum. 

A. Generation of parameter points 

The package SARAH was used to create a SPheno executable specific to the CBLSSM, which 
was used with the SSP package [75j to perform random scans over two different parameter 
regions. The first scan, which we refer to as the "democratic" scan, took random values for 
each diagonal entry of the R-sneutrino ^ bilepton Yukawa coupling Yx independently over 
its range. The other, which we refer to as the "hierarchical" scan, kept the (1, 1) and (2,2) 
entries as 10~^ and 10~^ respectively. The hierarchical scan was motivated by the hierarchy 
of the quark and charged lepton Yukawa couplings, but is not favored for any other reason. 
The ranges of the CBLSSM parameters for each scan are shown in Tab. [Tij 

In general, picking a set of GUT-scale parameters is very unlikely to result in a potential 
which has a phenomenologicaUy acceptable minimum. Instead of wasting time searching 
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Parameter 


Common to both 


Ml/2 


100 
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Mo 


100 - 


- 3000 


^0 


-3000 - 


- 3000 


tan /3 


3 


-45 


mz' 


1500 


- 3000 


tan/3' 


1.0 - 


- 1.5 


Parameter 


Democratic 


Hierarchical 


yll 

X 


0.05 - 0.6 


fixed 10-3 




0.05 - 0.6 


fixed 10-2 




0.05 - 0.6 


0.1 - 0.6 



TABLE 11. Ranges of parameters used in generating the samples. All dimensionful quantities are 
to be read as in units of GeV. The signs of fi and ji' were fixed to both be positive. The democratic 
scan consisted of 3158 points, the hierarchical of 2302 points. 

the parameter space for such points, the strategy adopted in SPheno (and others, such 
as ISA JET [76J, SOFTSUSY [77] and SUSPECT [78j) is to instead fix some parameters (in 
our case, /x, S^, /x', 5^/) at the SUSY scale, defined as Qsusy = ^/^^^~J^'} by demanding 
that the tadpole equations are satisfied for VEVs given as input (indirectly, by specifying 
m^, tan /3, m^/, tan (3') This is sufficient to ensure that the chosen point is an extremum, but 
does not guarantee that the extremum is not a saddle point or maximum. Any points that 
were found to be saddle points or maxima (by the presence of tachyons) were discarded. 

B. Finding all the tree-level extrema 

As mentioned in the introduction, finding the global minimum of even the tree-level 
potential as a function of several scalar fields is highly non-trivial. The Grobner basis 
method can be used to find all the minima of a polynomial function, as has recently seen 
some discussion in the literature 180] . A less well-known method is that of homotopy 
continuation [HH |82] , which has found use in several areas of physics [83l - [85] , in particular 
finding string theory vacua ^86j and extrema of extended Higgs sectors ^87j, where the 
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authors investigated a system of two Higgs doublets with up to five singlet scalars in a 
general tree-level potential. In contrast, the Grobner basis method is deemed prohibitively 
computationally expensive for systems involving more than five or six degrees of freedom 
[79] . while for our purposes eight to ten degrees of freedom are necessary, and the solutions 
for thousands of parameter points needed to be calculated. 

The numerical polyhedral homotopy continuation method is a powerful way to find all the 
roots of a system of polynomial equations quickly [88j . Essentially it works by continuously 
deforming a simple system of polynomial equations with known roots, with as many roots as 
the classical Bezout bound of the system that is to be solved {i.e. the maximum number of 
roots it could have) . The simple system with known roots is continuously deformed into the 
target system, with the position of the roots updated with each step. We made use of the 
publically-available program H0M4PS2 |[88j to find all the extrema of the tree-level potential 
via the polynomial tree-level tadpole equations. 

1. Vacuum expectation value phase rotations 

Degenerate lines of minima, corresponding to unconstrained phases of some of the VEVs, 

cause problems since the homotopy continuation method maps discrete numbers of solutions 

of a simple system to discrete numbers of solutions of the target system. From the physics 

point of view this is equivalent to the statement that only phase differences are relevant but 

not the phase of a single parameter itself. Already at the level of the superpotential in eq. [T] 

it is easy to see that two VEVs in the extended Higgs sector can be chosen real without 

changing the phase of a single parameter. We chose to take the VEVs of Hu and fj as purely 

real, and all other phases to be relative to thes^ As we are neglecting the fiavour structure 

of Yjy^ we can also choose all VR^i to be real (or equivalently we could take all VL^i real). 

If this fiavour structure is taken into account, only one of the six sneutrino VEVs can be 

taken real in general. However, even assuming that we can neglect fiavour mixing in y^, we 

have five real and five complex VEVs, resulting in fifteen real parameters which have to be 

calculated by H0M4PS2. Even though it is relatively fast, solving a system with ten degrees 

of freedom in typically twenty minutes on a single 3.4 GHz Intel 17 core, it takes about 10 

^ In the three generation case one could in principle also take in addition the VEVs of of and r] to be 

real be redefining the phases of /i and /i^ However, as this cannot be done consistently within H0M4PS2 

we take them to be complex and keep the two parameters real. 
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days for the 15 degrees of freedom. For this reason we decided that for the democratic scan, 
aU the VEVs would be taken to be real. Six points out of the 3158 were found which had 
one or more mass-squared that was negative at the deepest extremum found, and these are 
taken to have deeper minima where at least one of the VEVs is complex. In all six cases at 
least one sneutrino VEV was non-zero. The hierarchical case can effectively be treated like 
a one generation model and here we allowed for complex VEVs. We found in this case that 
every minimum with a non-zero imaginary part for any VEV was degenerate with another 
minimum of the same parameter point with purely real VEVs, with the same magnitudes but 
some different signs. Hence our conclusions would be unchanged if we restricted ourselves 
to purely real VEVs in the hierarchical scan. 

C. Loop corrections 

Before reaching any conclusion we have to estimate the nature of loop corrections to the 
scalar potential. The homotopy continuation method allows us to find all the solutions to 
the tree-level minimization conditions. However the effective potential at the one-loop level 
already has a more complicated structure including logarithmic contributions. 

Loop corrections might modify our previous discussion by changing the relative depth be- 
tween i?-parity-conserving and i?-parity- violating minima or promoting a gauge conserving 
minimum to be the global minimum at loop level. It might also happen (although less likely) 
that that the nature of the extremum is changed, e.g. a tree-level saddle point becomes a 
lower-lying minimum of the effective potential. 

In order to assess the effect of loop corrections, we evaluated the one-loop effective po- 
tential at all the extrema found by H0M4PS2 for every parameter point. In its more general 
form, the expression for the effective potential contains field-dependent masses. As our study 
centers around minima for the potential, the field-dependent masses are just the tree-level 
masses as function of the VEVs. We chose to remain in the Landau gauge and DR' scheme 
and follow the results in [89j . 

To evaluate the effective potential for a given set of VEVs we used the expression 

i ^ 

where is the tree-level scalar potential and i goes over all the real scalar, Weyl fermion 
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and vector degrees of freedom, with di accounting for all the spin, symmetry, and color 
factors. The rrii are the tree-level masses for a given set of parameters and VEVs at the 
renormalization scale Q. 

When the potential is non-convex the conventional loop expansion becomes complex [90]. 
Specifically, the value of the potential on a saddle point can be interpreted as an amplitude 
with a corresponding decay width (its complex part) [91j. It was sufficient to use the real 
part of the effective potential to estimate whether a saddle point might have been promoted 
to a lower-lying minimum by comparing it with the value of the effective potential in all of 
the other minima. 

By doing this it was possible to determine whether the hierarchy for the minima found at 
tree level was respected once loop corrections are included. We found that loop corrections 
have a non-negligible impact to the tree-level conclusions. 

IV. RESULTS 

We categorize the parameter points of our scans as in Tab. [TTT| Every point fell into one 
of these three categories: e.g. there were no points which had a global minimum with no 
unbroken symmetries, and there were no points which broke i?-parity without breaking both 
SU{2)l and U{1)b-l- 

As it happens, there were no "RPC" parameter points in either scan at either tree or 
one-loop level which had global minima deeper than the input minima, by which we mean 
the set of VEVs chosen as input for SPheno. 

The majority of "RPV" points broke i?-parity through non-zero sneutrino VEVs, but 
there were seven points in the hierarchical scan which had one-loop global minima with 
negative stop mass-squared, hence we assume that such parameter points have non-zero 
stop VEVs. We plot them along with the other "RPV" points since they do break i?-parity, 
in addition to SU{3)c and [/(l)e^. 

The third category, "unbroken", only appears in the one-loop-level analysis. It is known 
that particularly precarious arrangements of parameters that break symmetries sponta- 
neously at tree level may not break these symmetries when evaluated at one-loop level [92j . 
In such cases, the (not obvious) combination of Lagrangian parameters that parameterizes 
t^(l)5-L-breaking for example may be of the same magnitude as the loop corrections. We 
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Categorization 


Description 


"RPC" 


SU{2)l, U{1)b-l both broken, i?-parity conserved. 


"RPV" 


SU{2)l, U{1)b-l both broken, i?-parity broken. 


"unbroken" 


Either SU{2)l or U{1)b-l broken but not both, i?-parity conserved. 



TABLE III. Categorization of parameter points according to the symmetries broken by their global 
minima. 



Categorization 


Hierarchical scan 


Democratic scan 


total 


2302 


3158 




tree level 


one-loop level 


tree level 


one-loop level 


"RPC" 


1981 


1039 


3008 


2754 


"RPV" 


321 


358 


150 


131 


"unbroken" 





898 





267 



TABLE IV. Number of parameter points in the various categories. All of the parameter points from 
the hierarchical scan categorized as "unbroken" broke SU{2)l without breaking U{1)b-l im- 
parity, while 250 of the 267 parameter points from the democratic scan did so, with the remaining 
17 breaking U{1)b-l without breaking SU{2)l. Not all parameter points that are "RPC" at the 
one-loop level were "RPC" at tree level, and likewise for the "RPV" category. 

note that all such points that we found had zero sneutrino VEVs. 
A. Comparison with previous work 

Previous studies on how much of the BLSSM parameter space conserves i?-parity [3l] 
considered the positivity of all the soft SUSY-breaking mass-squareds of the R-sneutrinos to 
be the necessary and sufficient condition. However, as can be seen in Figs. [2] and |4| this is 
neither necessary nor sufficient. There are, of course, obvious trends, and parameter points 
with very low or negative m?c are as likely to break i?-parity as conserve it, but there are 
clearly parameter points which conserve i?-parity despite having a negative m?c. This can 
be understood by the fact that one has to consider the eigenvalues of the mass matrix which 
can be all positive because of contributions of the F-terms, D-terms and other soft SUSY- 
breaking terms for non-zero bilepton VEVs despite negative m~c entries. For instance, the 
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FIG. 1. Projections into various parameter planes of the 2302 hierarchical scan parameter points, 



categorized by the nature of their global minima (see Tab. Ill) at tree level. "RPC" points are 
plotted in green (light grey), and "RPV" points are plotted in red (medium grey). In the plots 
on the left, the "RPC" points are plotted on top of the "RPV" points, which are faded (very 
light grey), while in the plots on the right the "RPV" points are plotted on top of the "RPC" 
points, which are faded (very light grey). (m?c is the lowest or most negative of the three soft 
SUSY-breaking mass-squared parameters for the R-sneutrinos, evaluated at the SUSY scale.) 
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FIG. 2. Projections into various parameter planes of the 2302 hierarchical scan parameter points, 
categorized by the nature of their global minima (see Tab. |III[ ) at one-loop level. The scheme is as 
m Fig. [ll with now "unbroken" points plotted in blue (dark grey) and on top of the "RFC" points 
on the right, and faded and beneath the "RFC" points on the left. 
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FIG. 3. Projections into various parameter planes of the 3158 democratic scan parameter points, 



categorized by the nature of their global minima (see Tab. Ill) at one-loop level. The scheme is 
the same as in Fig. [Tj 

(z>3, element of the tree-level scalar mass-squared matrix is given by 



m. 



+^ {99BL -^l- ^l) + oIl (2 {vl - vl) + 3vl - f i) ) . 



(10) 
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FIG. 4. Projections into various parameter planes of the 3158 democratic scan parameter points, 



categorized by the nature of their global minima (see Tab. Ill) at one-loop level. The scheme is 
the same as in Fig. [2j 



Parameter points that have all positive m?c yet still have i?-parity- violating global minima 
are less surprising, as the presence of trilinear terms obviously could dominate for reason- 
ably large values of sneutrino and bilepton VEVs. Indeed, non-zero Vd and Vu are even 
occasionally necessary for non-zero bilepton VEVs in certain regions of parameter space. 
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where U{1)b-l is broken by bilepton VEVs despite the bilepton submatrix of the scalar 
mass-squared matrix having no negative eigenvalues for Vd = = 0. 

Furthermore, [34] concluded that a hierarchical Yx would always break i?-parity, because 
of the way the couplings enter the RGEs and drive the soft SUSY-breaking sneutrino mass- 
squared parameters negative before the bilepton mass-squared parameters, unless all three 
generations contribute substantially. The large fraction of "RFC" points in our hierarchical 
scan refutes this claim. Such points are harder to find (as evidenced by the smaller size 
of the hierarchical scan to the democratic scan), and in some sense rather finely-tuned, as 
evidenced by the large fraction of "unbroken" points where loop corrections to the diflFerences 
in potential are as large as the differences leading to the breaking of the symmetries at tree 
level. Still, almost half of the hierarchical points conserve i?-parity even at the one-loop 
level. 

V. CONCLUSION 

We have investigated the stability of i?-parity in a phenomenologically interesting region 
of the parameter space of the BLSSM by considering all the extrema of the tree-level potential 
and using the full one-loop evaluation of the potential at its extrema near those of the tree- 
level potential. We found that a significant fraction of points that were chosen to conserve 
i?-parity while having the correct mass for the weak vector bosons actually have global 
minima which violate i?-parity through non-zero sneutrino VEVs. However, we also found 
that more parameter points conserve i?-parity than violate it. 

We also found that both the conservation and violation of i?-parity were possible all 
over the parameter space, in contradiction to previous studies, and while trends in the 
parameters leading either case are visible, it is difl&cult to identify in advance what values 
of which parameters will lead definitely to the conservation or violation of i?-parity. 
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